Zero-Conductance Resonances due to Flux States in Nanographite Ribbon Junctions 
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The electronic transport properties through junctions in nanographite ribbons are investigated 
using the Landauer approach. In the low-energy regime ribbons with zigzag boundary have a single 
conducting channel of edge states. The conductance as a function of the chemical potential shows a 
rich structure with sharp dips of zero conductance. Each zero-conductance resonance is connected 
with a resonant state which can be interpreted as the superposition of two degenerate flux states 
with Kekule-like current patterns. These zero-conductance dips are connected with a pronounced 
negative magneto resistance. 
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The discovery of fullerene molecules and carbon 
nanotubes has triggered intensive research on various 
nanometer size carbon materials In these systems, 

the geometry of sp 2 carbon networks crucially affects the 
electronic states near the Fermi level. Studies with scan- 
ning tunneling microscopy and spectroscopy have con- 
firmed the connection between the electronic states of 
the single wall carbon nanotubes (SWCN) and their ge- 
ometry Recently the electrical transport measure- 
ment of individual SWCN became possible 0-0] and the 
quantized conductance of multi-wall carbon nanotubes 
was observed This initiated theoretical studies de- 
voted to effects of non-magnetic impurities Jf| , electron 
correlation and topological defects |pd|-|l5|. Besides 
these closed carbon molecules, there are also systems 
with open boundaries which display unusual features con- 
nected with their shape. They include small scale sys- 
tems based on graphite, so-called nanographites. There 
arc two basic shapes of regular graphite edges, called 
armchair and zigzag edges, depending on the cutting di- 
rection of the graphite sheet. Properties originating from 
such edges have been studied recently using the model of 
graphite ribbons, one-dimensional graphite stripes of in- 
finite length and finite width. Ribbons with zigzag edges 
(zigzag ribbons) possess electron states localized near the 
edge with energies very close to the Fermi level [l8| . 
Such states are absent for ribbons with armchair edges. 
The edge states of zigzag ribbons were analyzed in terms 
of nearest-neighbor tight binding models |l6|-[l8| and den- 
sity functional approach |19|. It was also pointed out 
that the edge states play important roles in the magnetic 
properties in nanometer-size systems due to their rela- 
tively large contribution to the density of states at the 
Fermi energy |l6|,[l8| . Recently, experimental evidence for 
edge states has been reported for nanographite systems 
derived from graphitized diamond nanoparticles |2(J ■ 

In this letter, we investigate the transport properties 
of ribbons related to these edge states. For this purpose 
we design three different types of junctions which con- 
nect zigzag ribbons of the same or different widths. The 
electronic states are described by a single-orbital nearest- 



neighbor tight binding model. The conductance of the 
junctions is evaluated using the multi-channel Landauer 
formula 1231, 



G(S) = ^^|V(S)| I 



(1) 



where t^ v (E) is a transmission coefficient from fi-th chan- 
nel to iv-th channel at energy E, calculated by a recursive 
Green's function method |24j]. Before discussing the de- 
sign of the junctions and their conductance properties, we 
review a few facts concerning the low-energy states in two 
kinds of graphite ribbons, the zigzag and the "bearded" 
ribbon (a variant of a zigzag ribbon). 

Zigzag ribbon: The zigzag ribbons are metallic for ar- 
bitrary ribbon width. The most remarkable feature is 
the presence of a partly flat band at the Fermi level, 
where the electrons are strongly localized near the zigzag 
edge. Each edge state has a non-vanishing amplitude 
only on one of the two sublattices, i.e. non-bonding 
character. However, in a zigzag ribbon of finite width, 
two edge states coming from both sides, have finite over- 
lap. Because they are located on different sublattices, 
they mix into a bonding and anti-bonding configuration. 
In this way the partly flat bands acquire a dispersion 
[Fig. 1(b)]. Note that the overlap is increasing as k de- 
viates from 7r/a, because the penetration depth of the 
edge states increases and diverges at k = 27r/3a, where 
a is the lattice constant. The dispersion depends on the 
ribbon width N (number of zigzag lines from one side 
to the other). Close to k — n/a, the spectrum has the 
approximate form = ±2tND^~ 1 (l — D]~/2), where 
Dk = 2cos(-y l ), and t is the hopping matrix element. 
Thus, although the edge states on each side separately 
have non-bonding character, together through their over- 
lap they provide one conducting channel except at ex- 
actly E = 0. The energy region of single-channel trans- 
port is restricted by the energy gap (A^) to the next 
channel A z ~ 4tcos[(7V - 1)tt/(2N +1)] @. 

Bearded ribbon: The bearded ribbon has one zigzag 
edge and one edge which has additional bonds (beard) 
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attached to the zigzag edge [Fig. 1(a) region M]. The edge 
states of both sides reside on the same sublattice so that 
in spite of their overlap their non-bonding character is 
retained and a completely flat band at E = is resulting 
for any width TV [Fig. 1(b)]. The absence of dispersion 
leads to the insulating behavior for the edge state chan- 
nel. The gap (Ah) to the first conducting channel is given 
by A 6 ~4icos[JV7r/(2/V + 2)] f|. 

We now turn to the design of a junction connecting 
two zigzag ribbons denoted by L (left) and R (right). As 
a first example we consider the case depicted in Fig. 1(a) 
where the junction region denoted by M is a bearded 
ribbon of length I (number of attached bonds). This 
junction model represents a metal-insulator-metal junc- 



tion providing an illustrative example for the peculiar 
transport properties of graphite ribbons regardless of the 
question as to whether bearded ribbons could be realized 
in nature. We calculate the conductance G(E) numeri- 
cally within the energy range \E\ < A&/2. The result for 
N = 30 and different values of I is shown in Fig. 1(c) using 
a logarithmic scale on the energy axis. The most striking 
feature is the large number of zero-conductance dips at 
energy values which depend on I. These dips represent 
resonances of total reflection. In the whole energy range 
no transmission resonances are observed for any value of 
I. Note also that the change of N does not lead to a 
qualitative change and only modifies the energy range of 
single-channel conductance. 
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FIG. 1. (a)The structure of the graphite ribbon junction. The length of the junction, /, is given by the number of attached 
bonds in the shaded region, (b) The schematic figure of energy band dispersion near E — of zigzag ribbons (left) and bearded 
ribbons (right), (c) The energy dependence of the conductance, for I — 5, 10, 15 and N = 30. (d) The energy dependence of 
the response functions, Xv an d Xv f° r I ~ 15 and N = 30. (e) The 3D plot of the correlation of circular current pattern in 
the M-region. (f) The energy dependence of the phase of the transmission coefficient of I — 15 when N = 30. 



Resonances are associated with discrete quantum lev- 
els in the junction regions. In the following we would 
like to characterize the quantum states encountered here 
and discuss the origin of total reflection. It was no- 
ticed earlier that some graphite ribbons form a triangular 
Kekule pattern of circular current driven by special ex- 
ternal boundary conditions Jl8[ . To investigate whether 
the junction states responsible for the resonance display 
a similar structure, we introduce the current vortex am- 
plitude defined on the dual (triangular) lattice as the 
clockwise circular sum of the currents Ii^ p on the bonds 



of the p-ih hexagonal plaquette, V p = S<=i ^i,p- An inci- 
dent current J m coming from lead L yields a current flow 
in the junction region, whose circulating component (vor- 
ticity), we probe by the following two response functions 

Xv = (l^pD/^in and x[? — {Vp)/Jim where the bracket 
(A) denotes the average value of A in the junction re- 
gion. The first function is a measure for the presence of 
circular currents and the second for the overall vorticity 
of the system. In Fig. 1(d) the energy dependence of both 
functions is shown for I = 15 and N = 30. The suscepti- 



2 



bilities Xv' 2 ^ show sharp divergences in the sense of linear 
response to the driving current J ln when E approaches 
the energy values of a zero-conductance dip. (Note that 
Xy' 2 "* are not defined at zero-conductance energies, be- 
cause no driving current J ln can flow.) This suggests that 
the states at these energy values form current vortex pat- 

(2) 

terns. The overall vorticity measured by Xv Jm changes 
sign at each resonance. We find that the circular current 
form a Kekule-like pattern. This is most easily seen in the 

Fourier transform F(k) = J2 P (j S ~J e 1 ^ Vp in the junc- 
tion region which has a clear peak structure of the trian- 
gle correlation of the circular currents pattern (Fig. 1(e)) 
(r p : coordinate of the ring center; k x (k y ): wave num- 
ber along (perpendicular to) zigzag lines in the junction). 
This type of pattern appears in F(k) very close to every 
zero-conductance dip and gradually disappears when de- 
viating from the resonance energies. Because F(k) is the 
correlation among the circular currents on the hexago- 
nal plaquettes, the state driven by the incident current 
is in this sense a flux state. The quantum state associ- 
ated with the resonance is the superposition of the two 
flux states with opposite circular currents, which com- 
bine to a standing wave within the junction region under 
the condition of time reversal symmetry. 

The presence of this quantum level plays an impor- 
tant role for the realization of the zero-conductance res- 
onances. At each resonance the transmission is not only 
carried by the usual tunneling through the insulating 
junction region, but also through this resonant state. 
This decomposition into two channels yields the follow- 
ing form for the transmission amplitude between the two 
leads close to a resonance at energy E , 



t(E) = t(E) 1 - 



iT/2 



E-E +iT/2J E-Ey 



t{E){E-E ) 



•tT/2' 
(2) 



where t(E) is a regular complex function of E and T 
the width of the resonance. The destructive interference 
of the two channels is a consequence of the symmetric 
form of the S'-matrix which relates the in- and outgoing 
waves of the two leads with and of the junction region 
p5| . This situation has been discussed in the context 
of the three-way splitter and a rigorous proof was given 
for the exact cancelling at the resonant energies p5j. A 
consequence of this resonant form of the transmission 
amplitude (Eq.(2)) is that the phase 8 of t(E) exhibits a 
jump by 7r at each resonance (see Fig. If) (similar to Ref. 
p6| , P7| ). The behavior of the phase is well-described by 
Eq.(2). The form of the S'-matrix mentioned above as 
well as the degeneracy of the flux states of opposite chi- 
rality, combining to a standing wave state, are based on 
time reversal symmetry. This symmetry can be destroyed 
by applying a magnetic field. We observe a pronounced 
negative magnetoresistance at the zero-conductance dip 
where the conductance grows proportionally to B 2 for 
small magnetic fields B, as we will discuss in detail else- 



where. We would like to mention here that circular cur- 
rents associated with zero-conductance dips were recently 
also reported for two-dimensional wave guides including 
junctions with stub geometry ]2q] . Also in this case a 
negative magnetoresistance should occur. 
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FIG. 2. (a) The structure of graphite ribbons junction type 
A and type B. (b) The energy dependence of the conductance 
for Nl = 50 and N R = 30. 

Next, we discuss two junctions with a more realistic 
design which do not contain bearded ribbons. They are 
shown in Fig. 2 (a) and connect zigzag ribbon leads of dif- 
ferent width (we show Nl — 50 and Nr = 30 as a rep- 
resentative case). The M-region contains a tilted zigzag 
edge for junction A and an armchair edge for junction B. 
The conductance of the two junctions is qualitatively dif- 
ferent. The tilted edge in A supports an edge state, sim- 
ilar to the bearded ribbon, on the same sublattice as the 
edge state on the other side. Indeed we find a very similar 
behavior of the conductance with large number of zero- 
conductance resonances [Fig. 2(b)] which are associated 
with flux states too. For junction B there is no localized 
state near the armchair edge. The conductance is rather 
featureless without any resonance in the energy regime 
of single-channel transport. There is a zero-conductance 
dip above the single-channel regime of the lead on the 
left-hand side (E > A Z (N = 50)/2 = 0.062t) for both 
junctions. The origin is also connected with analogous 
junction states which are based on flux states as the 

(1 2) 

analysis of the susceptibilities Xv ' show, although the 
conditions are different here, since the transmission oc- 
curs from three channels on the left hand side to a single 
channel on the right hand side. We omit here a detailed 
discussion of this more complicated case. For energies 
above the single channel threshold of the right-hand side 
lead (E > A Z (N = 30)/2 = O.lOli) no zero-dip features 
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appear m any case. 

From this analysis we conclude that the edge structure 
of the junction plays an important role in forming the 
resonant states in the M-region. It is well known that 
the electron states display chiral properties in graphite 
sheets, if there is an imbalance between the two sublat- 
tices e.g., by different onsite potential [p9| . The edge 
states on zigzag edges are a consequence of this imbal- 
ance, since the outer most sites belong to a single sublat- 
tice. The abrupt change of the sublattice on the edge as 
it occurs for previous junction and junction A yields the 
boundary condition to form the degenerate flux states 
in the M-region. These flux states combine to one res- 
onant state of standing-wave nature. Further numerical 
analysis shows that even single non-magnetic impurities 
and other simple defect structures disturbing the sublat- 
tice balance can cause a zero-conductance dip associated 
with flux states in zigzag ribbons. 

In conclusion, we numerically analyzed three types of 
nanographite junctions. We found that the conductance 
of various junctions having zigzag edges shows many zero- 
conductance dips as a function of energy (chemical po- 
tential). These dips are identified as resonances con- 
nected with resonant states based on flux states which 
form circular-current Kekule patterns. It is obvious that 
the topology of the edges is crucial for this phenomenon 
and the chirality connected with the sublattice structure 
plays an important role. The time reversal symmetry 
necessary for the existence of these resonances is violated 
by external magnetic fields, leading to negative magneto 
resistance. While the structures used in the calculation 
might be difficult to produce at present our results 
also suggest that transport properties of defective car- 
bon nanotubes, carpet-roll or papier-mache structures 
p2| could be rather different from the transport prop- 
erties of usual closed multi-wall nanotubes or SWCN 
which have only weak features in the low-energy regime 
Jl4| . The present study not only clarifies the importance 
of the edges and their shapes on transport properties, 
but also indicates the importance of theoretical studies 
to explicate the interplay between the transport prop- 
erties and the network topology of carbon atoms. The 
present numerical work provides the foresight concerning 
the analysis based on a low-energy effective theory, which 
is beyond the scope of this paper and will be presented 
elsewhere. Such an effective theory will serve as a basis 
for designing carbon-based electronic devices and for fur- 
ther theoretical work on effects of impurities or electron 
correlation. 
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